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Calculations have been made concerning the momentum that can be given to a jet of 
water by mixing it with a jet of high-speed steam. Apart from an application of the equa- 
tions of momentum and enthalpy, this raises questions concerning the speed of condensation 
and the acceleration of the water. 



1. Introduction 

Basic equations for the change of temperature 
and of momentum to be expected when jets of steam 
and water are mixed have been given by Thiamin 
and Dressier [l]. 3 

The present authors extended these calculations 
to a case where momentum must he communicated 
to a jet of water by mixing it with a jet of high- 
speed steam in such a way thai complete condensa- 
tion of the steam will take place. It was thought 
that condensation under constant pressure would have 
certain advantages. On the one hand there is a 
much smaller loss of useful momentum, while on the 
other hand tin 1 thermodynamic aspect of the con- 
densation becomes simpler. It was also thought 
that there would be less chance for the appearance 
of shocks. In order to obtain condensation under 
constant pressure, the chamber in which the condensa- 
tion takes place must have a sectional area decreasing 
in the direction of the (low (as is applied in the 
injector for steam locomotives), according to some 
law that must be determined. 

It was thought useful to have a concrete example 
in mind with definite numerical values. It was 
assumed therefore that the pressure in the mixing 
chamber should be kept equal to 0.8 atm. 

The water is supposed to enter into the mixing 
chamber in the axial direction, with a velocity 
V — lOm/sec and a temperature of T' = 20° C. The 
enthalpy of the water then is H' = 20.l kcal/kg, 
the density p' = 998 kg/m 3 . 

Data have been taken from W. Koch and K. 
Schmidt, VDI-Wasserdampftafem (Springer- Verlag, 
1952); since these data are presented in technical 
units, we must make use to a certain extent of the 
system in which the kilogram is taken as the unit 
of weight or force. We shall write kg* in this case, 
and reserve the notation kg when it refers to mass. 

The steam used as a driving force is supposed to 
he obtained from a boiler where the temperature is 
about 105° C and the pressure 7.5 atm (saturation 
pressure for this temperature). After coming from 
the boiler the steam is superheated to T = 360° C, 
the pressure remaining p =7.5 atm. The steam 
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then has an enthalpy H = 759A kcal/kg. It will 
he assumed that the expansion of the steam from 
7.5 atm to 0.8 atm occurs without friction losses 
nix! without change of entropy. If losses and 
entropy change cannol he avoided, we suppose 
tluit they are compensated for by a suitable in- 
crease of T and j) , so that we shall still have the 
same data for the expanded state as are given 
below. Aiter expansion the steam will then have 
the following properties: 

pressure p"=0.8 atm=8,000 kg* m 2 
temperature V v/ =99.0° C 
density p" = 0.462 kg/m 3 
enthalpy Z7"=639.4 keel kg 
velocity tf" = l,003 m sec. 

It is assumed that also the velocity u" of the steam 
is directed axially. 

The entrance area of the mixing chamber will be 
denoted by A; of this the part A s is taken in by (be 
water jet, and the part A (Is) by the steam jet. 

We write A m for the sectional area of the mixing 
chamber at the point where complete condensation 
has been reached; u m for the speed of the water 
when all is condensed; p m for the density of this 
water, for which we shall use a mean value Pw =985 
kg in 1 ; further T m for its temperature, H m for its 
enthalpy. We note that the mechanical equivalent 
of heat in technical units is J=427 kg*m per kcal; 
and that g= acceleration of gravity = 9.8 m/sec 2 . 

2. Basic Equations for the Mixing Process 

The following equations hold for the mixing- 
process : 
Equation oj continuity: 



Asp'u'+A (1—s) P " u"=A m Pm u n 



(1) 



When numbers are inserted, we solve for (A m /A) u m 
and obtain: 

(A m /A) w,„ = 0.47 + 9.56 s, (2) 

Equation of momentum: 

As P 'u n +A (1-s) P " u" 2 =A m Pm uj. (3) 
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Inserting numerical values we find: 

(AJA) u m 2 =4:72-372 s. (4) 

Equation of enthalpy: 

As P 'u'(H'+u' 2 l2gJ)+A(l-s)p"u"(H"+u" 2 l2gJ) 

=A m Pm u m (H m +u m 2 /2gJ). (5) 

Inserting numerical values we arrive at : 

(AJA) u m (H m +u m 2 /2gJ) = 25S-15d s. (6) 

From eqs (4) and (2) we can obtain u m ; from (6) and 
(2) we obtain H m -\-u m 2 /2gJ and then II m itself; from 
the latter quantity the temperature of the water can 
be found. It is also possible to find (AJA). 

The condensed water has the speed u m and the 
pressure p m —p' = 8,000 kg*/m 2 . When it is to be 
used for propulsion, it may have to be driven out 
against a higher pressure p e . In such a case it will 
be appropriate after the mixing to make the water 
flow through an expanding duct, in which the velocity 
decreases from the value u m to a value u e , adjusted so 
that the pressure will increase from p m to p e . In 
what follows we have assumed: P e = 2.0 atm= 
20,000 kg*/m 2 . We then find: 



Ue 



-u m 2 -239 (m/sec) 2 . 



(7) 



Formulas for other values of p e can easily be derived. 
The equations can be worked out for various 
values of s. After some trials it was found that 
s=0.5 gave a convenient case. Larger values of s 
(meaning that more water is admitted in the mixing 
chamber in comparison with the steam) make the 
exit speed rather low; smaller values of s (meaning 
that more steam is used in comparison with the 
water) make the final temperature higher and may 
lead to incomplete condensation. With s=0.5 the 
following results are obtained : 

u m =54.5 m/sec 
u e =52.3 m/sec 
H m =52.9 kcal/kg 
T m =52.9°C 

AjA=o.om. 

To the temperature T m corresponds a saturation 
pressure of 1,450 kg*/m 2 . Since the pressure in the 
mixing chamber is 8,000 kg*/m 2 , it was thought that 
there would be no danger for incomplete condensa- 
tion provided the time for condensation is not made 
too short (see section 4). 

3. Momentum and Efficiency of the Issuing 

Jet 

The total mass flow of water in the issuing jet is: 

A m p m u m kg/sec =A m y m ujg technical units of 

mass per second. (8) 



Here y m is used to represent the specific weight of the 
water, which, when expressed in kg*/m 3 , is given by 
the same number 985 as for its density in kg/m 3 . 
When the jet is applied to drive an underwater 
object with a speed of V m/sec, the usual arrange- 
ment will be that the water is taken out of the sur- 
roundings, where it had the speed V relative to the 
object, while it is reissued as a jet with the speed u e 
relative to this object. Hence the thrust to be ob- 
tained is given by: 

thrust = (A m y m ujg) (u e - V) kg*. (9) 

The useful work done in unit of time is: 

power= (A m y m ujg) (u e —V) Fkg*m/sec. (10) 

It may be possible to choose V so that this expres- 
sion has a maximal value; this occurs when V= J u e . 
We then find : 

thrust = \ (A m y m ujg) u e kg* (9a) 

power=i (A m y m ujg) u\ kg*m/sec. (10a) 

Inserting numerical values with 6—0.5, we arrive at: 

thrust = approx 13,800 A kg* (9b) 

power=approx 360,000 A kg*m/sec. (10b) 

In these formulas A is expressed in m 2 . When the 
entrance area of the mixing chamber is taken as 
100 cm 2 =0.01 m 2 , the values become: 

thrust = approx 138 kg* (9c) 

power = approx 3,600 kg*m/sec. (10c) 

In order to obtain a rough idea concerning the 
efficiency of the arrangement, we will consider only 
the energy used for the steam production. We have 
assumed that the water needed for this purpose is 
taken from the surroundings, at the temperature of 
20° C, so that the enthalpy is 20.1 kcal/kg. It 
follows that per kg water to be turned into steam, 
the enthalpy difference H -H' = 759A— 20.1 = 739.3 
kcal is needed. Since the mass of the steam needed 
per second is given by A (1 — s) p" u"=A (1— s)X463 
kg/sec, the energy required per second is : 

A (l-s)X463X 739.3X427 kg*m/sec, 

With ,9=0.5 this becomes: ^4X0.73 XI 8 kg*m/sec, 
and for ^4=0.01 m 2 the energy required (at least) is 
730,000 kg*m/sec. This would give a maximal 
efficiency of 

3,600/730,000 = 0.0049. 

The low value of the efficiency 7 should not be used 
as a measure for the merits or demerits of the system, 
since its merit is in the simplicity of operation. 
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4. Speed of Condensation 

For complete condensation not only a sufficient 
difference must be maintained between the pressure 
in the mixing chamber and the saturation pressure 
at the temperature after the condensation, but there 
is involved also a question of time. 

It is desirable that the condensation should require 
not too much time, since otherwise a mixing chamber 
of considerable dimensions will be needed in view of 
the velocity of the water which rises from u'=10 
m/sec to u m =54:.5 m/sec. In order to have a definite 
case it was thought that each element of volume of 
the water should stay not longer than 0.1 sec in the 
mixing chamber. The problem is whether complete 
condensation will be possible under this condition. 
Because we have condensation of steam in the 
presence of water, the question of obtaining a suffi- 
cient number of condensation nuclei does not present 
itself, but the important mutter is to have a sufficient 
surface area of the water in order to make rapid 
condensation possible. 

The speed of condensation depends on two factors: 

(a) Whether the difference between the pressure in 
the steam and the saturation pressure at the water 
surface will accelerate the steam sufficiently rapidly 
in order that it can reach the water surface within 
the time limit set above; and (b) whether the heat 
liberated by the condensation can flow off towards 
the interior of the water sufficiently rapidly. The 
latter problem depends on heat conduction and a 
simple estimate gave the impression that this was 
the really important question. 

In treating the problem of heat conduction into 
the interior a mathematical difficulty presents itself 
in consequence of the circumstance that the con- 
densation adds water to the already existing surface, 
so that the surface where condensation occurs dis- 
places itself with time. For the ease of a plane 
surface a treatment will be given (a) based on the 
ordinary solution of the one-dimensional heat equa- 
tion, in which this displacement is neglected; and 

(b) on a solution in which the displacement is taken 
into account according to a method indicated by 
I. I. Kolodner [2]. 

The equation for the heat conduction in one di- 
mension is, 






d 2 r 

*y 2 



(ii) 



where y is the distance measured from the surface 
of the water inwards. Here K=\/pC, where X is 
the thermal conductivity of the water; p is the 
density and C is the specific heat. For X we can take 
(approximately) 1.4X10 -3 cal sec" 1 cm" 1 °C~ l ; in 
c-g-s units p=l and C=l, so that k has the value 
1.4 X10" 3 cm 2 seer 1 . 

When the initial temperature of the water is 
T'(=20° C) and the surface on which the steam 
condenses acquires a temperature T s (=93° C, the 
saturation temperature for the pressure of 0.8 atm), 
the appropriate solution of the equation is, 



T=T'+a f 1-erf -^=Y (12) 

When we neglect the displacement of the surface, we 
must have T=T S for y = 0, which gives, 



a=T,-r / = 73°C. 



(13) 



The heat flow to the interior then becomes (per cm 2 
and per sec) 

2= - x (i)o=' C(T «- r,) V^ (14) 

The heat flowing to the interior in the time / is, 

2 — _ 

q.t=-r pC (T-D Vk*=*3.08V* cal/cm 2 . (15) 

Taking £=0.1 sec we find, 

gr |= 0.975 cal/cm 2 . (15a) 

The amount of heat to be taken up by the water 
is found as follows: There enters into the mixing 
chamber per second, A(l-s) p"u" kg steam. 4 

This steam is condensed into water with a temper- 
ature of at most 53° C (for s = 0.5). We take as 
mean temperature of the water: (20°+53°)/2=36° .5. 
The steam lias the enthalpy #" = 639.4 kcal/kg. 
The heat to be taken up by the water entering into 
the mixing chamber is given by the loss in enthalpy 
plus kinetic energy of the steam on condensation. 
Since the enthalpy plus kinetic energy of a given mass 
of steam remains the same as its value in the boiler, 
say 759.4 kcal/kg, and since the kinetic energy of the 
water is 0.1 kcal/kg, the rate at which heat is taken 
up by the water is 

A(\-s)p"u n (759.4-36.5-0.1) 

=4X0.5X0.462X1,003X722.8 

^1.675X10 8 A kcal/sec = 1.675 XI 8 A cal/sec. 

If we divide this by <Z we find that per second a sur- 
face of 

1.675X10 8 X^4/0.975X1.72X10 8 ^ cm- 
must be available. 

There enters the mass of water, 
A s p' u' kg/sec-5,000 A kg/see = 5X10 6 A r//see. 
which has the volume 5X10 6 

A cm 3 /sec (p= 1,000 kg/m 3 and w'=10 m/sec). 

It thus follows that per cm 3 of water we must have 
a surface area of 

1.72X10 8 /5X10 6 =34.4 cm 2 /cm 3 . 

When the water is dispersed into drops of radius R, 
the ratio of area to volume is 3/i?. Hence we re- 
quire drops of a radius not larger than 0.104 cm, 
that is, drops of a diameter of about 2 mm. 

4 In this formulas A is expressed as m 2 . However, this does not effect the re- 
sults since A finally drops out. 
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When the displacement oj the surface must be taken 
into account, we assume (following Kolodner [2]) 
that the surface at time t finds itself at a position, 



X=2b^Kt. 



(16) 



We now must have T=T S for y =—X=—2b -\fd. 
Instead of (13) this gives, 



T s -T' = a(l+erib). 



(i7: 



A second condition relating a and b is found by con- 
sidering the balance between the heat liberated by 
the condensation and the heat flowing to the interior; 
it can be written, 

where £=#"-36.5 = 603 kcal/kg=603 cal/g. 

With (12) and (16) we obtain 



\a J\ 

v^ exp( - 62)=pi6 V" 



The time factor drops out, and with \ = p Ck we find, 



a=^be b2 ~- 



Substituting this into (17) we arrive at 

C{T-T') 



&g»-(l+erf &)=- 



L\TT 



(19) 



(20) 



Xow with the data used above 

(7(T s -r)/ZV^=0.068; 
thus 

6 = 0.064, and a = 68.8° C. 

From (18) we now obtain 

g t =2 P LC-yj r id^2.90^t cal/cm 2 . 

Again taking £=0.1 we find 

^=0.916 cal/cm 2 , 

instead of 0.975 cal/cm 2 as given by (15a). 

It appears from these calculations that the effect 
of the condensation of water upon the surface makes 
the heat flow to the interior somewhat smaller, 
although the difference remains below 10 percent. 

Since it is necessary to disperse the water into fine 
drops in order to guarantee a sufficient condensation 
area, the calculations must be repeated for a spher- 
ical surface. This has been done only for the case 
where the displacement of the surface is neglected, 
in view of the mathematical difficulties involved in 
a more complete calculation. In this case it appears 



that with spherical drops the heat flow to the in- 
terior is smaller than with a plane surface. Probably 
the effect of the condensation upon the surface will 
again give a further reduction by at most 10 percent. 
Although the calculations cannot easily be made 
precise, they make it possible to obtain an estimate 
for the total surface area needed for the water 
entering into the mixing chamber per second, in order 
to allow the condensation of all the steam entering 
in the same period, when it is assumed that each 
water droplet will stay in the mixing chamber not 
over 0.1 sec. Fine dispersion of the water is neces- 
sary, in such a way that the drop size should be of 
about 2 mm in diameter. When this cannot be 
obtained, much more time must be allowed for the 
condensation. 

5. Equations of Motion for the Water Drop- 
lets 

Provisional equations have been developed for the 
motion of the water droplets in the mixing chamber, 
in order to find their speed as a function of the 
distance traversed. This is needed for the determi- 
nation of the cross-sectional area which the mixing 
chamber must have at various distances from its 
entrance. 

For this purpose it is assumed that the water has 
already been dispersed into fine drops and it is 
supposed that condensation takes place on these 
drops, so that no new drops are formed. 

The speed of the water drops will be denoted by 
u' \ this speed increases with time, and we can also 
consider u' to be a function of the coordinate x 
measured along the axis of the mixing chamber. We 
then can write for the drops, dx=u' dt, and express 
their acceleration as, du'/dt=u' (du'ldx). 

The speed of the steam will be denoted by u" . 
Also this speed can be a function of jc. 

We assume that the mixture of steam and water 
retains a constant pressure, and that the steam 
retains a constant temperature T" . The density of 
the steam, p" , will then be constant. The tempera- 
ture of the water droplets, on the other hand, 
increases in consequence of the condensation of 
steam upon them. The density of the water drops 
p will be considered as constant. 

We shall write A x for the sectional area of the 
mixing chamber at the position x. Let the number 
of water drops per unit volume at xben; the number 
of drops crossing a section of the chamber per second 
is then given by, 

N=A x nu'. (21) 

As mentioned it is supposed that no new drops are 
formed; when the state of motion is stationary the 
quantity N must be a constant (independent of t 
and x). 

The volume of a drop will be written, w=4:TrR s /3. 

The mass of steam condensing on a single drop in 
unit time will be denoted by M, which may be a 
function of x. We then have, 



p 'dw/dt=M. 



(22) 
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The speed of the drops increases in consequence 
of the momentum gained from the condensing steam. 
This momentum is given by M u" per second. It 
is possible that there is also a frictiona] force / 
acting on the drops in surplus of the direct gain of 
momentum (it may be, however, that a closer 
analysis of the phenomenon of condensation will 
lead to the result that/=0, since the steam in the 
neighborhood of the drop will flow towards the drop, 
so that there will be no appreciable shearing effect 
at the surface of the drop). 

Keeping both the direct transfer of momentum 
M u" and the frictional force / in the equations, we 
obtain for a single drop, 

p' d(wu')ldt=Mu"+f, (23) 

from which, with the aid of (22), 

P ' w du'/dt=M (u"-u') + f. (24) 

To form the conservation equations for the mix- 
ture of steam and drops we observe that the part of 
any section A x taken in by water is given byA z nv), 
while the part taken in by steam is A x {\— nw). 
Complete condensation will have occurred when nw 
has become equal to unity. 

The equation of continuity takes the form, 

. 1 riwp'ti'+A x (l-nw)p"u" = c(n\Mi\ni = K l . (25) 

7 

The equation of momentum becomes (since the 
pressure is constant), 

A x nwp'u' 2 +A x (l-nw)p"u'' 2 =conat&nt=K 2 . (26) 

The equation of enthalpy transport can be added to 
these equations; it can be us(mI to find the tempera- 
ture of the drops. We shall omit it here. 
From eqs (25) and (16) we deduce, 

4^pV(/-w')=if/-Z 2 . (27) 

With the aid of (21 ) this can be written, 

N P 'w(u"-u')=K 1 u"-K 2 . 
Differentiation with respect to the time gives, 

A 7 P '(u" '-u')(dw/df)+N p f w(du" '/dt)-N P 'w(du' '/dt) 

=K l {du"ldt), 

from which, 



{K l -Np f w){du ff ldt)=-Nj. 



(28) 



It follows that the speed of the steam u" changes 
only if there is an extra frictional force/. 



Equations (25), (26), and (21) can be used to 
express A f and n as functions of the other variables, 

x ~ P 'p"u'u"(u"-u') [ ' 

Np'p"u"(;u"~u') 



' (p'u' - P "u")K 2 -(p'u' 2 -p"u" 2 )K l 

(K l u"-K 2 )p"u" 

~'w[{p'u'- P "u")K 2 -{p'u' 2 -p"u" 2 )K y 



(30) 



The central problem is situated in the following 
three equations, 

dwldt = M/p' 

du'ldt \\l{u"-u') + 1\lp'w 

du"ldt=-j/(K 1 IN-p'w). 

If we attempt to neglect /, the last equation gives 
m" = constant. We then obtain from the other two 
equations, 

dwjdu' = w\ (u" — u') , 

the integral of which is, 

w(u" — u' ) = const an t . 

If we know the volume of the individual water 
drops at the entrance of the mixing chamber (a;=0), 
and if we make a scale for u' (we know that u f 
increases from the value 10 m/sec to the value 54.5 
m/sec with s = 0.5) we can find the course of w as a 
function of u' . 

Equation (30) then gives the corresponding value 
of //. When nw=l } the process is completed, 
Finally eq (29) would give the cross-sectional area 
A x as a function of u\ 

What remains is to find the relation between u' 
and x. This is determined by the equation 

du'/dx=(du'/dt)/u'=M(u''-u')/p'u'w, 

where du f /dt is obtained from (23) with /=0. In 
order to integrate this equation it is necessary to 
know M, and this must follow the analysis of the 
condensation process. 
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